Two-Dimensional Quantum Ferromagnets 
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We present 1/N Schwinger boson and quantum Monte Carlo calculations of the magnetization 
and NMR relaxation rate for the two-dimensional ferromagnetic Heisenberg model representing a 
quantum Hall system at filling factor v = 1. Comparing the analytic and numerical calculations, we 
find that the SU(JV) version of Schwinger boson theory gives accurate results for the magnetization 
at low temperatures, whereas the O(N) model works well at higher temperatures. 
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Two-dimensional quantum magnets have received par- 
ticular attention in recent years because of advances in 
materials synthesis associated with high-temperature su- 
perconductivity, thin films and surfaces, and semicon- 
ductor quantum wells. It has recently come to be appre- 
ciated that two-dimensional electron gases in quantum 
wells subjected to strong magnetic fields in the quantum 
Hall regime are novel itinerant ferromagnets. The strong 
external magnetic field quenches the orbital kinetic en- 
ergy but couples only very weakly to the spin degrees of 
freedom allowing low energy spin fluctuations to survive. 

Two-dimensional ferromagnets exhibit novel topologi- 
cal defects referred to as skyrmions by analogy with the 
corresponding objects in the Skyrme model in nuclear 
physics. What is unique about quantum Hall ferromag- 
nets is that these defects carry fermion charge and 
hence their ground state density can be controlled by 
moving the filling factor away from v = 1. The combina- 
tion of low energy spin fluctuations and these topological 
defects dramatically alters the NMR spectrum |3|,|| and 
the specific heat j|. 

With the advent of NMR and magnetoabsorption mea- 
surements in quantum Hall systems, it is now possible 
to measure the temperature dependence of the electron 
magnetization. Recent theoretical work j|] has evalu- 
ated the magnetization of this quantum critical system at 
v = 1 using SU(iV) and 0(N) formulations of mean field 
theory (with N = oo). In this paper we present analytic 
results for the 1 /N corrections to the magnetization and 
compare them with extensive quantum Monte Carlo sim- 
ulations. We also present numerical results for the NMR 
relaxation rate 1/Ti. In addition to their relevance to 
2D quantum ferromagnets in general and QHE magnets 
in particular, these results provide new information on 
the level of accuracy of large N expansion methods and 
show a highly non-trivial difference in the behavior of the 
SU(AO and 0(N) models. 

Low energy spin fluctuations in quantum Hall ferro- 
magnets at v = 1 are expected to be well described by 
the Heisenberg model || . The large N approach to this 
model is a systematic expansion around a mean field the- 
ory for N = oo and has the advantage of being equally 



valid at all temperatures T. Furthermore, even at the 
mean field level this approach correctly captures the fact 
that arbitrarily small thermal fluctuations destroy the 
long range order in two dimensions. An alternative mi- 
croscopic approach which includes spin-wave corrections 
to the electronic self-energy has also recently been devel- 
oped Trumper et al. |?|] employ an SU(2) Schwinger 
boson theory to describe Gaussian fluctuations in a frus- 
trated 2D antiferromagnet. 

We start from the Heisenberg Hamiltonian 
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(ij) 



on a square lattice, where S-S = 5(5+1) at each site. By 
introducing two Schwinger bosons at each site, 5 + = a r 6, 
5~ = b^a, and 5 Z = (a^a — wb)/2, subject to the con- 
straint a^a + b^b = 25, one obtains an equivalent boson 
Hamiltonian. The SU(2) spin algebra of this model is 
generalized to SU(A) After going over to the contin- 
uum, the SU(iV) Hamiltonian is written in terms of N 
bosons b a per site, 



H= d 



js(d,bi)(d 3 b a ) - ^bUdjb^bpid^) 
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where summation over repeated indices is implied and a 
is the lattice constant. The bosons are subject to the 
constraint ^ Q b' a b a = NS, which picks out the allowed 
spin states for total spin 5. The matrix h% — S a p(—l) a+1 
describes the coupling of the bosons to the magnetic field. 

The partition function can be written as a coherent 
state functional integral over complex fields b a . De- 
coupling the quartic term by means of a Hubbard- 
Stratonovich transformation one finds 



Z = J D 2 b a DXDQj exp^-i^ ^ dr J d 2 rC[b; 
with the Lagrangian density 

£ = 4 b* a dob a + JS {d 3 b* a )(d 3 b a ) + NJQ Q 
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+ UQj b* a {d 3 b a ) ~ UQj {db* a )b a - ^hf 3 bp a 
+ Xb* a b a -NSX, (3) 

where A is an imaginary Lagrange multiplier field which 
implements the constraint, and Q is the real- valued 
Hubbard- Stratonovich field which acts as a gauge field. 

The local equivalence of the groups SU(2) and 0(3) al- 
lows one to also write down an 0(3) boson model equiv- 
alent to Eq. ([!]), which can be generalized to 0(A) ||. 
The continuum 0(A) Hamiltonian reads 



H = 



d 2 , 
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with the constraints J2 a ^Ifia = NS/3 and J2 a ^a^a = 0- 
The matrix h contains A/3 copies of the 0(3) gener- 
ator matrix ((0, i, 0), (— i, 0, 0), (0, 0, 0)) along the diag- 
onal. The second constraint means that creating two 
bosons of kind, say, N produces a state which is a linear 
combination of states without two N bosons added. This 
constraint restricts the Hilbert space by identifying cer- 
tain states with one another. In the functional integral 
it is enforced by a new complex Lagrange multiplier /x, 
which enters the Lagrangian as //* b a b a /2 + [ib* a b* a /2. 

To obtain mean field results, the auxiliary fields A, Q, 
and, for O(N), fx are set to their values at the saddle 
point. We denote mean field values by a subscript 0. Af- 
ter Fourier transformation, the boson fields can be inte- 
grated out. The values of the auxiliary fields are obtained 
by solving the saddle-point equations. Because of gauge 
invariance, Qo can be taken to vanish. For O(A) we find 
fio — 0. The mean field magnetization Mq is then found 
from the free energy and has been given by Read and 
Sachdev for both models ||. 

Another useful quantity is the NMR relaxation rate 
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where A is the hyperfine structure factor, which we as- 
sume to be isotropic and momentum independent, Af is 
the number of sites, and S^q, ui) is the dynamic struc- 
ture factor, which is related to the transverse suscep- 
tibility by S(q,uj) = Im x + ~(q,w)/(l - e"^). The 
mean field rate is plotted in Ref. ||, here, we give an- 
alytic results for completeness. For SU(A), \/T\ = 
An B (Ao + [3B/2)/(32ir/3 2 J 2 S 2 ) and for O(A), 1/Ti = 
A [n B (A + 0B) + n B (A )} /(16tt/3 2 J 2 S 2 ). 

We now turn to 1/N corrections to the magnetization. 
For the SU(A) case, fluctuations in the auxiliary fields 
about their saddle-point values are denoted by iAX = 
A — A and AQ = Q, respectively. As a short-hand we de- 
note any fluctuation mode by re = {AA(r, r),AQj (r, r)}. 



Following the approach outlined in Ref. || , the partition 
function is written as a functional integral over the re, 



Z = Dr e e 



-NS[r e ] 



(6) 



with the action S = Sq + S<nr + Sioop > where 

S = A -1 Tr In Gq 1 , 

1 ( h(i f 
Sdir = J^J dr J d 2 r (NJQ • Q - NSX) , 

Sloop = N- 1 Tr In (1 + G v e r e ) . 

Here, the trace sums over momenta, Matsubara frequen- 
cies, and boson flavors a. Gg(k, iu> n ) = {—ihuj n + 
JSk 2 a 2 — Bha/2 + a 2 Xo)~ 1 is the mean field boson Green 
function, and the ve are vertex factors describing the 
coupling of the bosons to the fluctuation re. The ac- 
tion can be written as a Taylor series in fluctuations re, 

S = Si^Lo( 71 ') _1< ^li i Ti x " -rg n , where the n = 1 term 
vanishes since we are expanding around a saddle point. 
iSo conspires with the re independent part of Sdir to form 
the mean field free energy (3Fq = A£>(°) . The other terms 
in the series can be found by expanding 6>i oop and Sdir- 
The RPA fluctuation propagator is the inverse of the ma- 
trix s< 2 ). 

The magnetization is expressed in terms of boson oc- 
cupation numbers, M = l/N J2 a (^a^a), and the ex- 
pectation values (b^ba) are obtained by inserting a con- 
stant source term ^2 a jabjJ) a into the Hamiltonian in 
Eq. (||). We apply this procedure to the partition func- 
tion of Eq. (|^) and expand the result in a Taylor series 
in the re. The resulting Gaussian integrals can be evalu- 
ated by contractions over re. This corresponds to writing 
down all allowed diagrams with one external j a vertex 
and any number of internal vertices, connecting all in- 
ternal vertices by RPA propagators. Loops without an 
external vertex and with one or two internal ones are for- 
bidden since the former vanish in an expansion around 
a saddle point and the latter are included in the RPA 
propagator. Fig. [j] shows the diagrams of order 1/N. 

Summation over the frequency of the vertical RPA 
propagator in both diagrams must be done carefully, tak- 
ing into account normal ordering of operators at equal 
times. This procedure makes the frequency sums unam- 
biguous and removes a spurious divergence. The details 
will be given elsewhere |]l0| . The magnetization can now 
be calculated numerically. The momentum integrals have 
a logarithmic UV divergence, which is regularized by a 
lattice cutoff. The 1/N contributions from both AA and 
AQj decrease the magnetization, as intuitively expected. 

For the O(A) model we additionally have to deal with 
fluctuations A/i in the second Lagrange multiplier field, 
which couple only to bj^b^ and b a b a . To order 1/N the 
only contributions come from the two 1 /N diagrams with 
the vertical propagator replaced by a Afi RPA prop- 
agator and the direction of the boson Green functions 
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changed accordingly. The result converges for large mo- 
menta and actually increases the magnetization: Mean 
field theory, which enforces the constraint ^ Q b]j^ a = 
only on average, underestimates the magnetization be- 
cause it contains unphysical contributions from spin mul- 
tiplets of lower total spin. Analytic results are further 
discussed below. 

In order to test the accuracy of the analytic results, 
we have carried out Quantum Monte Carlo simulations 
using the stochastic series expansion method , which 
is ideally suited for the present calculation since it does 
not introduce any systematic errors. Sufficiently large 
lattices can be studied so that finite size effects are negli- 
gible. The method is based on a Taylor expansion of the 
density matrix e~^ H . Writing H is terms of its one- and 
two-body terms, H — Y]^f. 1 Hj, the partition function 
can be written as [p"lf 

a n=0 S„ ' i=l 

where S n denotes a sequence of indices (h, h, l n ) with 
k 6 1, ...,M, and |a) = |S* , Sf , . . . , Sjy) is an eigenstate 
of all the operators Sf. The relative weight in this expres- 
sion can be made positive definite by adding a suitable 
constant to H . For a system of finite Af and (3 only se- 
quences of finite length contribute significantly and the 
limit n — > oo poses no problem (the average power (n) is 
given by \E\[3, where E is the total internal energy). 

We want to emphasize an important feature that 
makes the sampling particularly efficient: The external 
field is chosen in the x direction. This automatically 
causes the simulation to become grand-canonical and 
there are no longer any problems associated with a re- 
stricted winding number. If the transverse field is not 
too weak {B/J > 0.02), it causes the auto-correlation 
times of all calculated quantities to become very short, 
even though only purely local updates are used. Fur- 
thermore, it enables easy access to observables involving 
both diagonal and off-diagonal operators. Details of the 
implementation will be presented elsewhere p| . 

For a 4 x 4 system we have compared our QMC data 
with exact diagonalization results, and they agree to 
within statistical errors. Relative errors are typically of 
the order 10 -4 for all system sizes considered. For all 
the field strengths presented in this paper the results for 
16 x 16 and 32 x 32 sites agree to this precision (finite-size 
effects increase with decreasing B), and we here present 
magnetization results for the larger size. 

The q dependent imaginary time susceptibility is 

X + -(q,r) = ^(^(r,r)5 z (0,0)) e ^ r . 

r 

Once it has been calculated, the dynamic structure factor 
can be obtained by inverting the relation 



X+-(q,r) = i f duiS&uy-r" 

using the maximum-entropy technique Jig ]. The B field 
introduces an additional complication to the analytic 
continuation since it introduces a delta function peak in 
the spectral weight at u> = B arising from the q = 
response. The maximum entropy method cannot resolve 
separate delta functions in the spectrum and since we 
are studying relatively small fields this potentially causes 
problems for the calculation of 1/T±. We therefore found 
it useful to separate out the q = sector before continu- 
ing the average X) q X H ( C 1? T )- The results are, however, 
still sensitive to the somewhat sharp peaks in the spectral 
function near u> = B from momenta near q = 0; a prob- 
lem which becomes worse with increasing system size. 
We therefore only use 16 x 16 sites for the 1/T\ calcu- 
lation, which should be sufficient in the interesting tem- 
perature regime. We estimate the statistical errors using 
the bootstrap technique. As with any results obtained by 
the maximum entropy method, the error estimates have 
to be viewed with some caution. 

Fig. ^| shows the magnetization for spin S — 1/2 and 
magnetic field B — 0.1 J. The results are typical for 
all fields considered (0.02 < B/J < 0.32). The SU(AT) 
results agree with Monte Carlo data only at the lowest 
temperatures. In fact, at low temperatures the SU(AT) 
mean field magnetization is known to agree with the 
non-interacting magnon approximation for the Heisen- 
berg model up to exponentially small corrections. This 
is not the case for the O(N) model. There is a distinct 
crossover to moderate and high temperatures, where the 
O(N) 1/N results agree quite well with the data, whereas 
the SU(AT) 1/N correction eventually becomes much too 
large. The main differences between the systems studied 
above are that (i) the analytic calculations use a con- 
tinuum approximation, and (ii) the leading terms in the 
1/N expansion do not include skyrmion effects, whereas 
the Monte Carlo simulations in principle do. Both effects 
should lead to discrepancies at temperatures of order J 
and higher. We will further pursue the question of why 
the 0(N) model works better than the SU(iV) model at 
most temperatures in a subsequent paper p0| . 

Experimental data by Barrett et al. Q show a more 
rapid drop off at higher temperatures, probably due to 
disorder affecting the normalization of the data. In Fig. || 
we compare our results with recent magnetoabsorption 
measurements by Manfra et al. fl4| ]. The experiments 
agree well with our 0(N) 1/N and Monte Carlo results 
except at low temperatures. The calculations have been 
done for B — 0.32 J, which is the experimental field cor- 
rected for finite width of the quantum well [|14| . The 
approximation used by Kasner and MacDonald |6j leads 
to a much higher magnetization than obtained here. 

In summary, we find from our comparison of Monte 
Carlo and analytic approaches that the SU(N) model cor- 
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rectly captures the low-temperature magnetization and 
that the first 1/N correction to the O(N) model gives 
accurate results at higher temperatures. We also find 
that the NMR relaxation rate evaluated by analytic con- 
tinuation from the Monte Carlo data agrees rather well 
with the mean field result, as shown in Fig. ^. Both the 
1/N and Monte Carlo calculations were rather involved. 
Technical details will be presented elsewhere |l^,^2[ • 
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FIG. 2. Magnetization of a 2D quantum ferromagnet with 
magnetic field B = 0.1 J and spin S = 1/2. 
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FIG. 1. Diagrams of order 1/N for Solid lines de- 

note boson mean field Green functions, wriggly lines are RPA 
fluctuation propagators, and dashed lines are external j a in- 
sertions. 
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FIG. 3. Magnetization for B = 0.32 J compared with ex- 
periments from Ref. [14J. 
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FIG. 4. NMR relaxation rate 1/Xi for B/J = 0.1. 
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